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Abstract 

The CLEO III experiment has recently accumulated a large statistics sample of 4.73 • 10 6 T(35) 
decays. We present the first evidence for the production of the triplet Y(1Z?) states in the four- 
photon cascade, T(35) -> 1X b{2P), X b(2P) -> jT{lD), T(1D) -> wCbO-P), Xb(lP) -> 7^(15), 
followed by the T(15) annihilation to e + e~ or fi + fi~ . The signal has a significance of 9.7 standard 
deviations. The measured product branching ratio for these five decays, (3.3 ± 0.6 ± 0.5) • 10~ 5 , is 
consistent with the theoretical estimates. We see a 6.8 standard deviation signal for a state with 
a mass of 10162.2 ± 1.6 MeV/c 2 , consistent with the T(1Z>2) assignment. 

We also present improved measurements of the T(3S) — > 7r°7r°T(lS') branching ratio and the 
associated di-pion mass distribution. 



*Submitted to the 31 st International Conference on High Energy Physics, July 2002, Amsterdam 
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I. INTRODUCTION 



Long-lived bb states are especially well suited for testing QCD via lattice calculations JTJ] 
and effective theories of strong interactions, like potential models ||. The narrow triplet 
S states, T(15), T(2S) and T(3S), were discovered in 1977 in proton-nucleus collisions at 
Fermilab [[J. Later, they were better resolved and studied at various e + e~ storage rings. Six 
triplet P states, x fe (2Pj) and x&(1Pj) with J = 0, 1, 2, were discovered in radiative decays of 
the T(3S) and T(2S) in 1982 || and 1983 [|], respectively. There has been no observations 
of new narrow bb states since then despite the large number of such states predicted below 
the open flavor threshold (see Fig. [I]). 

In this paper, we present the first observation of the T(1D) states. Except for the yet 
unobserved T(2D) states, these are the only long-lived L = 2 mesons in nature. They are 
produced in a two-photon cascade starting from the T(3S) resonance: T(3S) — > 7Xb(2Pj), 
Xfc(2Pj) — > jT(1D). To suppress photon backgrounds from 7r°s, which are copiously pro- 
duced in gluonic annihilation of the bb states, we select events with two more subsequent 
photon transitions, T(1D) — > 7x&(lPj), x fe (lPj) — > jT(1S), followed by the T(IS') annihi- 
lation to either e + e~ or (see Fig. [I]). The product branching ratio for these five decays 
in sequence was predicted by Godfrey and Rosner [|7| to be 3.76 ■ 1CT 5 . 

The CLEO III experiment has recently accumulated 4.73 -10 6 T(3S*) decays, which makes 
the detection of such rare processes possible. The CLEO III sample constitutes roughly a ten- 
fold increase in T(3S') statistics compared to the CLEO II data set ||, and roughly a four- fold 
increase compared to the integrated CUSB data set |{J. Thanks to the good granularity and 
the large solid angle of the CLEO CsI(Tl) calorimeter, the CLEO III detection efficiency 
for these final states is about 4.5 times larger than in the CUSB detector. Even though the 
CLEO III calorimeter is essentially the same as in the CLEO II detector jTU|, the photon 
selection efficiency and detector resolution were improved in the endcaps and in part of the 
barrel calorimeter by the installation of a new, lower- mass, tracking system |TIj. Another 
change was the replacement of the time-of-flight system by a RICH detector in the barrel 
part. Finally, the calorimeter endcaps were restacked and moved farther away from the 
interaction point to accommodate a new, higher luminosity, interaction point optics. 



II. DATA SELECTION 

We select events with exactly four photons and two oppositely charged leptons. The 
leptons must have momenta of at least 3.75 GeV/c. We distinguish between electrons and 
muons by their energy deposition in the calorimeter. Electrons must have a high ratio 
of energy observed in the calorimeter to the momentum measured in the tracking system 
(E/p > 0.7). Muons are identified as minimum ionizing particles, and required to leave 
150 — 550 MeV of energy in the calorimeter. Stricter muon identification does not reduce 
background in the final sample, since all significant background sources contain muons. Each 
photon must have at least 60 MeV of energy. We also ignore all photons below 180 MeV 
in the calorimeter region closest to the beam, because of the spurious photons generated by 
beam-related backgrounds. The total momentum of all photons and leptons in each event 
must be balanced to within 300 MeV/c. The invariant mass of the two leptons must be 
consistent with the T(IS') mass within ±300 MeV/c 2 . Much better identification of the 
T(1S) resonance is obtained by measuring the mass of the system recoiling against the 
four photons. The average resolution of the recoil mass is 17 MeV/c 2 . Cuts on the recoil 
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FIG. 1: 66 mass levels as predicted by one of the potential models. The levels are denoted by 
their spectroscopic labels, n 2S+1 Lj, where n is the radial quantum number (n = 1, 2, . . .), S is the 
total quark spin {S = spin singlets, S = 1 spin triplets), L is the orbital angular momentum 
(L = S, P, D, . . .) and J is the total angular momentum of the state ( J = S + L). The four-photon 
transition sequence via the T(1D) states is shown. An alternative route in four-photon cascade 
via the T(25) state is also displayed. 




mass are described below. The mass resolution of the produced T{1D) state depends on the 
measurement of the energies of the two lowest energy photons in the event. Thus, we require 
that at least one of them is detected in the barrel part of the calorimeter, where the energy 
resolution is best. The energy resolution as determined by the fit to the T(35') — > ^Xb (2-Pj) 
photon lines is a Ej = (4.6 ± 0.2) MeV for £ 7 = 100 MeV 03 . 



III. T(35) -» 7r vr T(15) TRANSITIONS 

The most prominent signal among T(3S') — ■> 7777T(1S") decays are T(3S') — > 7r 7r°T(15') 
transitions. To identify these decays, we try all three combinations of four photons into 
two-photon pairs. We find the minimal 7r° mass deviation chi-squared among these combi- 
nations: xlo n o = min {[(M^ 7 - M^)/a M a] 2 + [(M^ - M^o)/^] 2 }. After the require- 
ment Xw°7i-° < 6, the deviation of the four-photon recoil mass from the T(IS') mass, 
AM = M recoi ; 47 — Mr(is), is found. The distribution of the ratio, AM/<r(AM), where 
cx(AM) is the expected recoil mass resolution, is plotted in Fig. |[ As can be seen from the 
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sidebands of this distribution, backgrounds from non-T(3S') — > T(IS') transitions (e.g. radia- 
tive Bhabha and fi— pair events) are small. We subtract this background by fitting the recoil 
mass distribution with the signal shape extracted from the Monte Carlo simulations. The fit 
is also shown in Fig. 0. We observe 737±28 events with an efficiency of 13.6%. The efficiency 
was determined from the Monte Carlo simulations. This sample contains a small fraction, 
0.9%, of photon cascade background that we estimate from the cross-efficiency determined 
with the Monte Carlo (below 1%) and the branching ratios measured from the data, as de- 
scribed below. After the background subtraction, we obtain the following result for the prod- 
uct branching ratio: B(T(3S) -> 7T 7r°T(15))-i3(T(15) -> l+r) = (5.67±0.22±0.35) • 10" 4 . 
The first error is statistical and the second is systematic. The systematic error is dominated 
by uncertainties in the Monte Carlo simulation of the selection efficiency. It was determined 
by varying the selection criteria. Here, B(T(1S) — ► means either B(T(1S) — > e + e~) or 
B(T(1S) — ► since these two branching ratios should be equal from lepton universal- 

ity. Using the world average value 0, B(T(1S) — > = (2.43 ± 0.06)%, we can unfold 

our result for £(T(3S) -> 7r 7r°T(lS)) = (2.33 ± 0.09 ± 0.16)%. This result is consistent 
with, but more precise, than previous measurements by CLEO II [Q, (2.03±0.28±0.19)%, 
and CUSB ||, (2.3 ± 0.4 ± 0.3)%. The efficiency-corrected invariant mass distribution of 
the two-pion system is shown in Fig. ^. As in the previous measurements, some dynamical 
structure is observed in this distribution. 




AlWa(AM) M(7t°7t ) (MeV) 



FIG. 2: Distribution of the recoil mass devi- 
ation from the T(1S) mass for the T(35) — > 
7r°7r°T(lS') candidate events. The solid-line 
represents the fit of a signal peak on top of 
the linear background (dashed-line). 



FIG. 3: Efficiency-corrected distribution of 
the 7r°7r° invariant mass observed in the data. 



The determination of B(T(3S) —> 7r 7r T(lS')) presented above serves as a good cross- 
check of our experimental techniques. In a separate paper submitted to this conference 
fl~2|| , we also present the analysis of the two-photon cascades between the T(3S f ) and the 
T(15) or T(2S). Branching ratios for two-photon cascades via x fe (2P) and x fe (lP) states 
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were determined with improved precision compared to previous experiments. These data, 
together with the 7r°7r° transitions, were used to obtain the calibration of absolute photon 
energies in the 60 — 800 MeV range. Using the well known masses of the T(1S), T(2S) 
and T(3S) together with the recoil mass distributions like the one shown in Fig. |2], we 
calibrated the photon energies to a precision of 0.35%. The photon energies observed in the 
two-photon cascade data for T(3S) — > jXb^Pj) agreed well with the previous measurements. 
Thus, the systematic effects in both rate and photon energy measurements are well under 
control. 



IV. EVIDENCE FOR THE T(LD) STATES 

To look for four-photon cascades between the T(3S) and T(1S), we must suppress back- 
grounds from the 7r°7r° transitions. We first require, X^o n o > 6. In addition, every photon 
pair must have an invariant mass at least two standard deviations away from the nominal 
7T° mass. 

To look for T(1P) events, we constrain events to be consistent with a photon cascade from 
the T(3S) to the T(1S) via one of the Xb(2Pj) and one of the Xbi^Pj) states (J = 0, 1,2). 



The masses of the P states are well known from other measurements For each J 2 p, J\p 
combination we calculate a chi-squared: 



xId,j 2P ,j 1p ( m t(id)) 




— _& 7J - {Mr {ID), J2P, JlP) 



cr Eyj 



where E 1 j are the measured photon energies; _g^ ected are the expected photon energies 
from the masses of the bb states and the measured photon directions in the cascade. The 
masses of the T(1D) states are not know. Therefore, we minimize the above chi-squared 
with respect to Mr (id) which is allowed to vary for each event. The above formalism 
requires that we know how to order the four photons in the cascade. While the highest 
energy photon must be due to the fourth transition, and the second highest energy photon 
must be due to the third transition, there is sometimes an ambiguity in the assignment 
of the two lowest energy photons to the first two transitions, since the range of photon 
energies in the T(35) — > 7Xb(2Pj) decay overlaps the possible photon energy range in the 
Xb{2Pj) — > 7T(1£)) transition. We choose the combination that minimizes the above chi- 
squared. There are nine possible combinations of J2P, J\p spins. We try all of them and 
choose the one that produces the smallest chi-squared, x\d = m ^ n x\D j 2P Ji P - 

In addition to the four-photon cascade via the T(1D) states, our data contain events 
with the four-photon cascade via the T(2S) state: T{3S) — > ■yXb(2Pj), Xb(2Pj) — > jT(2S), 
T(2S) -> iXb(lPj), Xb(lPj) -> lT(lS), T(1S) -> l+l- (see Fig. 0). This is a coin- 
cidence of the two-photon cascades, T(3S) —> ■yyT^S) and T(2S) — > jjT(1S), which 
were experimentally observed in the two-photon cascades detected in the T(3S) [|l^] and 
T(2S) data fp~5| . The product branching ratio for this entire decay sequence (including 
T(IS') — > l + l~) is predicted by Godfrey and Rosner [^] to be 3.84 • 10~ 5 , thus com- 
parable to the predicted T(1D) production rate. Combining the measured rates for 
B(T(3S) -> 77T(2S)) -B(T(2S) -> l+l~) @ and B(T(2S) -> jjT(1S)) ■ B(T(1S) -> l+l~) 
|15| with the measurement of B(T(2S) —> fl3fl , we obtain a factor 1.9 ± 0.3 higher 

value. In these events, the second highest energy photon is due to the second photon tran- 
sition (see Fig. p. Unfortunately, these events can be sometimes confused with the T(1P) 
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events due to the limited experimental energy resolution. The second and third photon tran- 
sitions in the T(2S) cascade sequence can be mistaken for the third and second transition in 
the T(1D) cascade sequence. Therefore, it is important to suppress the T(2S) cascades. We 
achieve this by finding the J2P, J\p combination that minimizes the associated chi-squared 
for the T(2S) hypothesis, x\s = m ^ n X2S,j 2P ,Jip^ wnere xls * s exactly analogous to Xid with 
the Mx(id) replaced with Mx(2S)- We then require xls > 12- Notice that the masses of all 
intermediate states are known for the T(2S) cascade, thus this variable is more constraining 
than Xid- To further suppress the T(2S) cascade events, we construct a quasi- chi-squared 
variable, xis^ that sums in quadrature only positive deviations of the measured photon 
energies from their expected values. This variable is less sensitive to fluctuations in the 
longitudinal and transverse energy leakage in photon showers that sometimes produce large 
negative energy deviations and correspondingly a large xis va hie. With the additional cut 
xls > 3, the cross-efficiency for T(2S) events is reduced to 0.7% and is further suppressed 
by a factor of two when a cut of x\d < 10 is required. 
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FIG. 4: xiv distribution for the data. The FIG. 5: xId distribution for the T(LD) 
cut value used for the mass analysis is indi- (solid-histogram) and T(2S) cascade Monte 



cated by the vertical bar. 



Carlo (dashed-histogram). The latter distri- 
bution was scaled up by a factor of 10 to be 
visible on the scale of the signal Monte Carlo. 



The data x\d distribution, after we require the four-photon recoil mass to be between —4 
and +3 standard deviations from the nominal Y(IS') mass, is shown in Fig. |j. A narrow peak 
near zero is observed, just as expected for Y(1.D) events. The signal Monte Carlo distribution 
for T(1D) events is shown in Fig. [5[ The background Monte Carlo distribution for the T(2S) 
cascades, after a factor of 10 enhancement relative to the T(1D) normalization, is also shown 
for comparison. Clearly, the T(2S) cascade background cannot produce as narrow a peak 
as observed in the data. The estimate of the number of T(2S) cascade events passing the 
Xid < 10 requirement is 1.6—3.0 events (depending on the assumed T(2S) product branching 
ratio). We observe 44 events after this cut in the data. The cross-efficiency for 7r°7r° events 
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is 0.02% in this sample, which gives a negligible contribution to the signal region. A flat 
tail in the Xid distribution for the data, not observed for the signal Monte Carlo, indicates 
that there is still some background present from other processes like radiative QED events. 
Since such processes can satisfy the multiple mass constraints built into the Xid variable 
only by a random coincidence, the distribution of that background varies slowly in \\d 
except for the bias towards lower values introduced by the \\d minimization with respect to 
Mx(id)i J2P and J\p. We use the T(3S) — > 7r 7r°T(lS') Monte Carlo (with the n° veto cuts 
removed) as a model of the Xid distribution for this background component (see Fig. ^|). 
Normalization of this background is fixed from the observed tail of the Xid distribution. 
This shape is likely to overestimate the amount of bias in \\ D towards the lower values 
since, unlike the T(35) — > 7r°7r°T(lS') process, the QED backgrounds do not automatically 
satisfy the Mt(3S) — Mr(is) mass constraint built into the chi-squared variable. Therefore, 
we also used a linear background model, with the line parameters fixed by a fit to the tail of 
the x\d distribution. Depending on the background shape model, the \id < 10 sample has a 
9.5 — 13.9% background rate (including the T(2S) cascade background). A fit of the T(1D) 
signal, on top of the backgrounds, is shown in Fig. |7|. The T(2S) background contribution 
was fixed from the Monte Carlo simulations. Normalization of the other backgrounds was 
left as a free parameter in the fit. The fit gives a signal amplitude of 40.7 ± 6.8 events. If 
we try to fit the data with the background contribution alone, the likelihood value of the 
fit drops dramatically. From the change of the likelihood we estimate that the signal has 
a significance of 9.7 standard deviations. Significant signals are observed in the 7777e + e~ 
and 7777/i + yU~ subsamples taken separately (5.2cx and 8.2<r, respectively). 

V. MASS ANALYSIS 

To turn the observed signal yield into a signal branching ratio, we need to estimate the 
signal efficiency Even though the signal shape in the \\d variable is fairly independent of 
the masses and spins of the produced T(LD) states, the selection efficiency does depend 
on these parameters via the cuts used to suppress the T(2S) cascade background. This is 
illustrated in Fig. |8|, where we plot the selection efficiency as a function of the mass of the 
T(1D) state. Three different curves are presented, for J 1D = 1,2 and 3. The J\o = 2 state 
is predicted by Godfrey and Rosner to be produced with the highest rate, 2.6 ■ 10~ 5 . The 
J2P = 1, J ip = 1 combination is predicted to be the main production and decay path for 
this state. The selection efficiency for the J\d = 2 state gradually increases with its mass, 
but is reasonably high in the entire search window. The J\d = 3 state, predicted with a 
0.8 • 10~ 5 rate, can only be produced via the J 2 p = 2, J 1P = 2 combination due to angular 
momentum conservation. This happens to be an unfortunate combination, since for a wide 
range of masses such a production fakes the T(2S) cascade, and the detection efficiency 
after the T(2S) suppression cuts is very low. We become sensitive to this state only if its 
mass exceeds about 10170 MeV/c 2 . Finally, the J\d = 1 state is predicted by Godfrey and 
Rosner to have the lowest product branching ratio, 0.4 • 10~ 5 . This state is again predicted 
to be produced mostly via the J 2 p = 1, J\p = 1 combination, and therefore our efficiency 
for this state is similar to the J\d = 2 efficiency 

A straightforward way to measure the mass of the produced Y(1.D) state is to calculate 
the mass of the system recoiling against the two lowest energy photons in the event. The 
distribution of the difference between the center-of-mass energy and this recoil mass, Ecu — 
M recoi i 11121 which should peak at My^s) ~ -^t(id), is shown for the Monte Carlo simulation 
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of a T(1D 2 ) state with a mass of 10160 MeV/c 2 in Fig. |9|. This distribution can be well 
described with a Gaussian that has a power-law tail on the low side (a so-called Crystal 
Ball Line Shape), as illustrated in Fig. |9|. The asymmetric tail is caused by transverse and 
longitudinal shower energy leakage in the calorimeter. The Gaussian part has a resolution 
of a= (6.8 ±0.3) MeV/c 2 . 

The mass value that minimizes xId * s 8 ^ so an estimate of the true T(1D) mass. The 
distribution of this mass for the same Monte Carlo sample as above is shown in Fig. [10]. A 
narrow peak at the correct mass is observed with a resolution of a = (3.1 ±0.1) MeV/c 2 . The 
mass resolution is improved in this method due to the constraints from the well measured 
masses of the P states built into the x\d variable. Unfortunately, there are also two satellite 
peaks, one on each side of the main peak. These false peaks originate from photon energy 
fluctuations which can make a wrong J 2 p, J\p combination produce the smallest chi-squared 
value. 

The distributions of E C m ~M recoi i 7l72 and the T(1D) mass which minimizes x\d f° r eacn 
event (called the most-likely mass) are shown for the data in Figs. |TT| - |rZ| after the x\d < 10 
cut. The background fraction for this sample is in the 9 — 14% range. Two peaks are 
observed in the distribution of the most likely mass, around 10160 and 10175 MeV/c 2 . The 
recoil mass distribution does not show a double-peak structure, but the mass resolution is 
worse in this variable, as discussed above. 

To analyze the most likely mass distribution, which has a better mass resolution, we fit the 
data to either a two-peak or a one-peak (positioned either around 10160 or 10175 MeV/c 2 ) 
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FIG. 8: Selection efficiency (with the Xid < 10 cu * included) for T(1D) states of different J\d 
as a function of the mass of the state. 



hypothesis. The background is assumed to have a flat distribution. Since the backgrounds 
are not significant, their parameterization is rather unimportant. Because, in principle, up 
to three triplet T{1D) states can contribute to our sample, we also tried to force a third 
peak into the fit function. However, the fit either makes the amplitude of this peak zero, 
with no effect on the amplitudes of the other two peaks, or the third peak collapses into one 
of the other two peaks. The two-peak fit, displayed in Fig. gives the best confidence level 
(58%). The amplitude (mass) of the first peak is 27.8^ events (10161.2 ± 0.7 MeV/c 2 ), 
and of the second peak 12.0^4;g events (10174.2 ± 1.3 MeV/c 2 ). The assumption that there 
is no peak around 10160 MeV/c 2 produces a very low confidence level (0.04%) and can be 
ruled out. From the change of likelihood between these two fits, the significance of the first 
peak is 6.8 standard deviations. The hypothesis of just one peak around 10160 MeV/c 2 
also produces a good confidence level of 43% (see Fig. [TJJ), with 38.6ig;2 events in the peak 
(10162.0 ± 0.5 MeV/c 2 mass). Furthermore, while the recoil mass distribution is consistent 
with the production of two states with the amplitudes and masses as determined by the fit 
above (confidence level 36%), it prefers to allocate most of the events to the lower mass state, 
37.7^5, with only 2. events attributed to the higher mass state, when the amplitudes 
are let free in the fit. This fit is displayed in Fig. [15] and has a confidence level of 48%. The 
single-peak interpretation of the recoil mass distribution also has a high confidence level 
(50%), with 40.5^| events at a mass of 10163.4 ± 1.3 MeV/c 2 (Fig. [HJ). 

In the region of the first peak in the most likely mass distribution, (77 ±9)% of the events 
have the J<ip = 1, J\p = 1 combination as the most likely. This is consistent with the Monte 
Carlo simulations for both the Y(1_D 2 ) (78%) and T(lDi) (73%) states. Furthermore, we 
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FIG. 9: Distribution of the difference be- 
tween the center-of-mass energy and the re- 
coil mass against the two lowest energy pho- 
tons for the T(1I?2) Monte Carlo events gen- 
erated with a mass of 10160 MeV/c 2 . The 
solid-line illustrates the fit of the Crystal Ball 
Line Shape to the distribution. 



FIG. 10: Distribution of the most likely 
mass assignment for the T(1Z?2) Monte 
Carlo events generated with a mass of 10160 
MeV/c 2 . In addition to the narrow peak at 
the correct mass, two satellite peaks at lower 
and higher mass are produced. 
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FIG. 11: Distribution of the difference be- 
tween the center-of-mass energy and the re- 
coil mass against the two lowest energy pho- 
tons for the data. 



FIG. 12: Distribution of the most likely 
mass assignment for the T(1D) candidate 
events in the data. 
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have very low efficiency for the T(1D^) state at this mass (Fig. [8]). When we remove the 
cuts suppressing the T(2S) cascades, which are responsible for the efficiency loss for the 
T(lDs) state, the amplitude of this peak does not increase dramatically. Therefore, this 
cannot be the Jm = 3 state. The Y(1_D 2 ) state is predicted to be produced with a 6 times 
larger rate than the T(1D 1 ) state The Clebsch-Gordan coefficients play a big role in 
this ratio, thus the model dependence of this prediction is weak. Therefore, we take the 
T(1D 2 ) state as the most likely interpretation of the first peak. Taking an average over the 
different fits described above, we measure mass of this state to be 10162.2 ± 1.6 MeV/c 2 . 
The error is dominated by the dependence on the fit approach (±1.2 MeV/c 2 ), with smaller 
contributions from the statistical error (±0.7 MeV/c 2 ), photon energy scale error (±0.7 
MeV/c 2 ) and the uncertainty on the T(3S') mass (±0.5 MeV/c 2 ). Predictions of the most 
successful potential models predict the center-of-gravity of the T(1D) triplet to be around 
this mass Calculations of the fine splitting of the T(1D) states predict the Jn> = 2 state 
to be close to the center-of-gravity mass 0. Our mass measurement is consistent with these 
predictions. 

Only 3 events have J 2 p = 2, Jip = 2 as the most likely combination. If the Jio = 3 
state was produced with the amplitude of 12 events, as suggested by the two-peak fit to the 
distribution of the most likely mass, we would expect 9.4 such events, including 1.5 events 
of feed-down from the Y(1D 2 ) state at 10162 MeV/c 2 . As described previously, the recoil 
mass distribution shows no evidence for a state at 10174 MeV/c 2 either. Thus, there is no 
strong evidence for the second state in our data. More data, which will be accumulated by 
CLEO III at the T(3S") in the future, will help to clarify the interpretation of these results. 

VI. INCLUSIVE PRODUCTION RATE FOR THE T(1D) STATES 

To calculate the signal efficiency in the fit to the Xid distribution (Fig. |7|), we weight the 
Y(1D 2 ) efficiency for a state with a mass of 10160 MeV/c 2 and the T(1D 3 ) efficiency for a 
state with a mass of 10175 MeV/c 2 in a ratio proportional to the observed peak amplitudes 
in the two-peak fit to the observed most likely mass distribution. We use the difference 
between this weighted efficiency and the individual efficiencies for each state as an estimate 
of the systematic error. This gives an efficiency of (13.2 ± 1.0)%. Converting the number 
of T(1D) events determined by the fit to the Xid distribution to a product branching ratio, 
we obtain (3.3 ±0.6 ±0.5) ■ 10~ 5 . This is an inclusive rate that sums over all possible T(1D) 
states. This rate is in good agreement with the 3.76 • 10 -5 prediction by Godfrey and Rosner 

VII. SUMMARY 

In summary, we present evidence at the 9.7 standard deviation significance for production 
of T(1D) states in four-photon cascades from the T(3S). The inclusively measured product 
branching ratio, B(T(3S) - jXb(2P))- B( Xb (2P) - 7 T(LD))- B{T{ID) -> 1X b(lP))- 
B(xb(lP) -> 7T(1S))- B(T(1S) -> l + l ) = (3.3±0.6±0.5) ■ 10~ 5 , is in good agreement with 
the theoretical predictions by Godfrey and Rosner J/J, where B(T(1S) — > above means 
the average branching ratio over e + e~ and fi + fi~ . From the observed mass distribution we 
see a 6.8 standard deviation signal for a state with a mass of 10162.2 ± 1.6 MeV/c 2 . This is 
likely to be the J\d = 2 state, though we cannot rule out the J\n = 1 hypothesis. However, 
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FIG. 13: The fit (solid-line) of a two-peak FIG. 14: The fit (solid-line) of one-peak 
structure to the observed distribution of the structure to the observed distribution of the 
most likely mass. The dashed-line indicates most likely mass, 
the fitted contribution of the first peak plus 
the background. The signal shape for each 
peak has a central peak, and a smaller satel- 
lite peak on each side of the central peak. 
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FIG. 15: The fit (solid-line) of a two- 
peak structure to the distribution of Eqm — 
Mrecou^^. The peak positions were fixed to 
the masses obtained by the fit to the most 
likely mass (the higher mass state appears 
at a lower value of Ecu — ^recoil 7172) • The 
amplitudes of the peaks were allowed to float. 
The dashed-lines indicate the fitted contribu- 
tions of each peak separately. 



FIG. 16: The fit (solid-line) of a one-peak 
structure to the observed distribution of the 
difference between the center-of-mass energy 
and the recoil mass against the two lowest en- 
ergy photons. The amplitude and peak posi- 
tions were free parameters in the fit. 
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the J\d = 1 hypothesis is strongly disfavored by theoretical expectations. Evidence for a 
second heavier state in our data is so far inconclusive. This is the first discovery of an L = 2 
quarkonium state below the open flavor threshold. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. We thank J. L. Rosner for the useful discussions con- 
cerning the T(1.D) states. M. Selen thanks the PFF program of the NSF and the Research 
Corporation, and A.H. Mahmood thanks the Texas Advanced Research Program. This work 
was supported by the National Science Foundation, and the U.S. Department of Energy. 



[1] See e.g. C. T. Davies et al. (UKQCD Collaboration), Phys. Rev. D58, 054505 (1998) and 
references therein. 

[2] For reviews see e.g. D. Besson and T. Skwarnicki, Ann. Rev. Nucl. Part. Sci. 43, 333 (1993); 

E. Eichten and C. Quigg, Phys. Rev. D49, 5845 (1994). 
[3] S. W. Herb et al, Phys. Rev. Lett. 39 (1977) 252; W. R. Innes et al, Phys. Rev. Lett. 39, 

1240 (1977) [Erratum- ibid. 39, 1640 (1977)]. 
[4] K. Han et al. (CUSB), Phys. Rev. Lett. 49, 1612 (1982); G. Eigen et al, Phys. Rev. Lett. 49, 

1616 (1982). 

[5] C. Klopfenstein et al, Phys. Rev. Lett. 51, 160 (1983). and Ref. |]. 
[6] F. Pauss et al (CUSB), Phys. Lett. B 130, 439 (1983). 

[7] S. Godfrey and J. L. Rosner, Phys. Rev. D64, 097501 (2001); see also W. Kwong and J. L. Ros- 
ner, Phys. Rev. D38, 279 (1988) and other references therein. 
[8] G. Crawford et al (CLEO II), Phys. Lett. 294B (1992) 139. 

[9] U. Heintz et al (CUSB), Phys. Rev. D46 1928 (1992). We corrected the B(T(3S) -» 
7r°7r°T(15)) published in this paper to the value of B(T(1S) — > l + l~) used in our calcula- 
tions. 

[10] Y. Kubota et al, Nucl. Instrum. Meth. A320, 66 (1992). 
[11] D. Peterson et al, Nucl. Instrum. Meth. A478, 142 (2002). 

[12] "Study of Two-Photon Transitions in CLEO III T(3S) Data", G. Bonvicini et al. (CLEO), 
paper submitted to this conference, ICHEP02 ABS949, CLEO CONF 02-07. 

[13] K. Hagiwara et al. (Particle Data Group) Phys. Rev. D66, 010001 (2002). 

[14] F. Butler et al. (CLEO II), Phys. Rev. D49, 40 (1994). We corrected the B(T(3S) -» 
7r 7r T(15')) published in this paper to the value of B(T(1S) — ► l + l~) used in our calcula- 
tions. 

[15] Ref. § and W. Walk et al. (Crystal Ball), Phys. Rev. D34, 2611 (1986). 



14 



